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ABSTRACT 

When  a nearly  fiat  fluid  surface  z = o is  subject  to  a small 
sinusoidal  distrubance  of  the  form  S'Z  z&z.9  cos  kx  and  when  gravity  - 

i or  an  equivalent  acceleration  field  - acts  on  the  fluid  in  such  a 

\ 

■,  sense  as  to  destabilize  the  surface,  then  the  irregularities  grow 

i 

exponentially  with  time  in  the  regime  of  small  amplitudes  where 
f Y linear  theory  applies  (Rayleigh).  At  larger  amplitudes  linear 

theory  does  not  apply  and  the  gas  or  magnetic  field  on  the  other 
side  of  the  interface  penetrates  linearly  with  time  into  the  fluid 
' in  the  form  of  fingers  or  bubbles  (G. I . Taylor).  Considered  here 

is  the  subsequent  stage  of  successive  mergers  between  these  fingers 
to  make  larger  and  faster  moving  fingers.  A primitive  and  tenta- 
tive analysis  of  this  merger  process  suggests  the  conclusion  that 
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the  advance  of  the  leading  - and  growing  - fingers  relative  to  the 
position  of  the  ideal  interface  is  described  by  an  acceleration 
which  is  of  the  order  of  1/50  of  the  acceleration  responsible  for 
the  instability.  Observations  are  not  available  to  check  this 
conclusion,  which  however  gives  some  indication  of  what  effects  to 
look  for  when  experiments  like  those  of  Lewis  and  Allred  and  Blount 
are  pushed  into  the  regime  of  very  high  accelerations.  To  the 
extent  that  the  statistical  analysis  is  valid,  it  suggests  that 
under  the  conditions  assumed,  a free  slab  of  fluid  cannot  be 
accelerated  to  more  than  roughly  50  times  its  own  thickness  with- 
out suffering  breakthrough. 

I.  INTRODUCTION  AND  SUMMARY 

Magnetic  fields  can  be  envisaged  today  with  a strength 
so  great  (much  more  than  10^  gauss)  that  the  resulting  magnetic 
pressure  can  accelerate  a metal  object  to  a high  velocity.  We  are 
concerned  here  with  a situation  where  the  field  is  so  intense 
(field  energy  per  unit  volume  of  the  order  of  heat  of  fusion  per 

zr 

unit  volume  or  B^lCrgauss)  that  the  metal  becomes  liquid  and 
the  conditions  are  at  hand  for  the  development  of  Rayleigh-Taylor 
instability  with  perturbations  in  the  interface  between  magnetic 
field  and  liquid  of  ever  Increasing  magnitude.  Under  the  extreme 
conditions  that  can  be  envisaged  no  way  presents  Itself  to  hold 
down  this  instability. 

The  surface  tension  of  metal  (Section  II)  is  high  but  the 
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accelerations  considered  are  so  extreme  (fa  ^ 10  \or  lO1^)  cm/sec2) 
that  surface  tension  prevents  only  the  growth  of  those  instabili- 
ties whose  reduced  wavelength  X = X/S.T  is  of  the  order  of  ^ /o~%rIc 
cm  or  less.  The  idea  is  discussed  eo  increase  the  effective  surface 
tension  by  covering  the  surface  with  an  alternating  laminar  struc- 
ture composed  <.f  two  fluids  of  high  mutual  surface  tension,  thus 
enhancing  the  effective  surface  tension  by  a factor  proportional 
to  the  number  of  lamina.  However  the  factor  of  stabilization 
achievable  in  this  way  is  not  great  enough  to  prevent  the  growth 
o^  interface  irregularities  with  reduced  wavelengths  of  micrometer 
and  greater  magnitudes. 

An  alternative  idea  is  discarded  in  the  present  context  with- 
out discussion:  to  reduce  the  rate  of  growth  of  interface  dis- 

turbance^/! th  X less  than  some  specified  value  Xt  ( by  inter- 
posing betv/een  the  principal  metal  and  the  magnetic  field  a layer 
with  thickness  ^ , in  which  the  density  falls  off  smoothly 
from  that  of  the  principal  metal  to  very  small  values^ \ We  limit 
attention  to  metal  sheets  so  thin  that  this  approach  is  not  feasible. 

Viscosity  is  a third  effect  which  acts  to  inhibit  the  growth 

of  Taylor  instability,  but  a simple  analysis  shows  that,  in  the 

extreme  circumstances  here  considered,  viscosity  is  ineffectual  in 

diminishing  the  rate  of  growth  of  instabilities  of  reduced  wave- 

il  -6 

length  X = 10“e%{or  IOcttJ  or  less. 

In  the  absence  of  any  effective  v/ay  to  inhibit  them  Taylor 
instabilities  will  grow  beyond  the  small  amplitude  regime  of 
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(l  2) 

exponential  rise'  * ' and  beyond  the  regime  of  formation  of 

spikes  and  fingers  - and  linear  rise  of  these  fingers  or  bubbles - 
considered  by  Taylor'  J and  will  come  into  a third  regime  where 
larger  bubbles  grow  by  capture  of  smaller  ones  (Fig.  l).  An  attempt 
is  made  in  Section  III  at  an  order  of  magnitude  statistical 
analysis  of  this  late  phase  of  Taylor  instability.  In  default 
of  relevant  experimental  information  about  the  details  of  this 
late  stage  the  theoretical  analysis  is  necessarily  uncertain.  In 
so  far  as  the  simplified  assumptions  and  approximations  used  in 
Section  III  are  valid,  they  would  suggest  that  the  head£  of  the 
bubbles,  or  fingers,  eat  their  way  through  the  liquid  with  an 
effective  acceleration  -relative  to  the  ideal  interface-  of  the 
order  of  1/50  of  the  bulk  acceleration  of  the  fluid  itself.  In 
other  words  under  the  conditions  assumed  it  is  concluded  that  it 
is  impossible  to  propel  a sheet  of  fluid  metal  to  more  than  rough- 
ly fifty  times  its  thickness  without  breakthrough.  This  estimate 
should  be  taken  not  as  a reliable  guide  for  the  design  of  a pro- 
pulsion system  but  as  an  indicator  for  the  kind  of  observations 
and  measurements  that  might  be  attempted  in  future  experiments. 

In  so  far  as  the  statistical  analysis  of  Tay3or  Instability 
has  any  validity  it  provides  a new  approach  to  an  old  question: 

How  far  does  it  pay  to  go  in  trying  to  produce  an  ideal  surface 
in  an  effort  to  reduce  Taylor  instability?  The  analysis  of  this 
point  (Section  IV)  suggests  that  it  is  impossible  ever  to  approach 
so  close  to  ideality  that  further  improvements  are  without 
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Fig.  1.  Later  or  "non-linear"  phase  of  Taylor  instability  (B  — >C) 
compared  and  contrasted  with  the  earlier  "linear"  phase  (A— HB).  In 
the  linear  phase,  the  shorter  is  the  reduced  wave  length  * of  a dis- 
turbance, the  greater  is  its  exponential  rate  of  rise:<*=  (g/^t )?. 
This  law  of  growth  ceases  to  be  valid,  however,  when  the  amplitude 
of  the  disturbance  has  become  comparable  with  % . Also  the  shape 
of  an  initially  exactly  sinusoidal  disturbance  of  the  surface 
ceases  to  be  sinusoidal  when  it  has  grown  into  the  non-linear 
regime . 
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significant  pay  off.  On  this  account  electro-polishing  is  briefly 
recalled  as  one  technique  by  which  to  improve  the  smoothness  of 
a surface. 

II.  THE  MINOR  EFFECTS  OF  SURFACE  TENSION  AND  VISCOSITY  IN 
INHIBITING  TAYLOR  INSTABILITY  AT  HIGH  ACCELERATION 

When  a fluid  of  density  and  surface  tension  o'  is 

plastered  smoothly  on  the  celling  of  a room  a small  disturbance 

in  the  surface 

S Z = SZ0  cert  Ax  (l) 

will  grow  exponentially  in  time  in  accordance  with  the  formula 

SZ  - «fZ00  €**  (2) 

or 

$Z  - SZ oo  . (3) 

In  a more  complete  discussion  one  allows  for  the  effect 
of  surface  tension  o'  (erg  / cm2): 

cC2  = gk  - ( cr /p  )r3  (4) 

or  the  effect  of  viscosity  (dynamic,  , g/cm  sec;  kinematic, 
z/=7 //,  cm2/sec)  or  both  in  inhibiting  the  growth  of  the 
disturbances  ^ ) . 

Experiments  have  confirmed  the  expected  effect  of  surface 
tension.  There  is  no  reason  to  doubt  that  the  expected  effect 
of  viscosity  will  appear  when  appropriate  experiments  are 
performed.  Both  effects  one  would  like  to  increase  dramatically 
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if  by  so  doing  the  growth  of  Taylor  instability  were  effectively 
inhibited. 


TABLE  I.  Reduced  wavelength  ^CRlT  above  which  disturbances  on 
water  and  copper  are  unstable  against  exponential  growth  of 
small  amplitude  disturbances.  From  Eq  (l4)  of  text,  K&Vt'“ 


Acceleration 

H2° 

Cu 

Surface  tension 
(g/secd) 

73  (20°C) 

1103  (1131°C) 

Acceleration  in 
Case  I 

10^  cm/sec2 

10^  cm/sec2 

^cRir , 1 

Acceleration  in 
Case  II 

2.7  x 10"^cm 
1013  cm/sec2 

3.6  x 1(5^  cm 
1013  cm/sec2 

2.7  x 10"6  cm 

3.6  x lO~^cm 

The  effective  surface  tension  is  non-isotroplc  when  the  push 
comes  from  a magnetic  field  with  lines  of  force  running  in  the 
y-directlon  parallel  to  the  surface  (yz-plane).  From  Faraday's 
picture  of  tensions  acting  along  the  lines  of  force  it  is  clear 
( Chandrasekhar ( 6) ) that  this  tension  augments  the  natural  surface 
tension. 
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whose  circular  wave  number  k = (o,  k&)  points  exclusively  along 
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the  magnetic  lines  of  force.  Under 'the  conditions  of  interest 
here  the  magnetic  fields  are  so  strong,  the  augmentation  of  surface 
tension  is  so  great  and  the  stabilization  thereby  achieved  is  so 
marked  that  surface  disturbances  running  in  this  direction  will 
be  disregarded.  On  the  other  hand,  when  the  disturbance  is 
described  by  a vector  k running  in  the  x-direction  the  magnetic 
lines  of  force  in  effect  behave  like  separate  bundles  of  rubber 
tubing,  every  other  bundle  being  raised  and  intermediate  bundles 
depressed  by  the  perturbation.  In  this  case  the  magnetic  field 
provides  no  assistance  in  stabilizing  the  surface  and  the  effective 
surface  tension  is  to  be  identified  with  the  surface  tension  of 
the  metal  alone.  To  this  extent  it  is  irrelevant  that  the  accelera- 
tion is  brought  about  by  a magnetic  field  instead  of  a gas.  In 
another  respect  the  difference  is  significant.  In  the  magnetic 
case  the  pattern  of  spikes  and  bubbles  will  b p drawn  out  along  the 
magnetic  lines  of  force  in  such  a way  as  to  create  not  spikes 
and  bubbles  but  crests  and  troughs.  However,  the  cross  section 
of  these  topologically  different  types  of  irregularities  we  envisage 
to  be  not  greatly  different.  Consequently,  the  factor  sv  1/So 
in  the  statistical  analysis  of  bubble  growth  is  thought  to  be 
changed  somewhat  in  absolute  value  but  not  in  general  Import  when 
one  goes  from  the  bubble  case  to  the  magnetic  case.  Therefore 
the  bubble  theory  developed  in  Section  III  will  not  be  revised 
on  this  account,  and  in  most  of  the  ensuing  discussion  the  aniso- 
tropy will  be  disregarded,  that  is  the  Taylor  instabilities  will 
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be  treated  as  If  they  arose  from  pressure  by  a gas.  This  is  pre- 
sumably the  case  in  which  experiments  can  be  carried  out,  to  judge 
from  the  work  of  Lewis ^ and  Allred  and  Blount where,  however, 
the  accelerations  were  too  small  and  the  surface  tension  too  great 
to  allow  opportunity  for  the  later  amalgamative  phase  of  T'^lor 
instability  to  develops.  In  contrast,  for  the  accelerations  listed 
in  Table  I the  natural  scale  of  dimensions  for  the  first  bubbles 
is  so  small  that  macroscopic  observations  will  reveal,  not  these 
first  bubbles,  but  only  the  much  larger  ones  that  arise  from  them 
by  many  generations  of  merger. 

In  principle,  one  can  greatly  increase  the  effective  surface 
tension  operative  at  t.ie  interface  between  two  fluids,  A (magnetic 
field)  and  C (molten  copper),  by  replacing  the  simple  interface 
by  a composite  interface  (Fig.  2)  built  out  of  many  thin  lamina 
of  liquid  C and  another  liquid  B (a  plastic,  for  example), 

according  to  the  pattern ... AAAACBCB CBCCCC.  Ordinarily  such 

a structure  is  difficult  to  produce  and  even  after  being  produced 
will  disappear  through  the  mechanism  of  formation  of  drops  and 
rearrangement  into  a structure  such,  for  example,  as 

AAAABBBCCCCCCC 

Fabrication  and  preservation  of  the  laminar  structure  would 
seem  much  simpler  in  the  case  of  a solid  metal  sheet  which  is 
destined  to  be  converted  into  a liquid  only  at  the  time  of  sudden 
buildup  of  a powerful  magnetic  field  just  above  its  surface.  Let 
the  homogeneous  metal  sheet  be  replaced  by  a metal  sheet  of  laminar 
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structure  produced  by  alternate  electroplating  in  two  electrolytes 
of  very  different  constitution  or  otherwise.  Let  the  two  substances 
be  so  chosen  that  their  interfacial  tension  is  as  large  as  possible. 
Then  the  effective  surface  tension  of  the  entire  molten  metal 
sheet  with  respect  to  deformations  of  reduced  wavelength  “k  of 
the  order  of  the  foil  thickness  or  greater  has  been  increased  by 
a factor  of  the  order  of  the  number  of  lamina. 

Consider  for  example  (Table  II)  a metal  sheet  of  thickness 
10“2  cm2  made  of  alternate  sheets  of  metal  and  plastic  each 
100  atoms  or.  ^ 10“^  cm  thick.  The  effective  surface 

tension  for  small  amplitude  disturbances  of  reduced  wavelength 

> 10"2  cm  will  be  increased  above  its  normal  value  by  a 

factor  of  the  order  ^ N = 10“  2 cm/lO“^  cm  = 10^;  thus 

•SSe&tefciVe Ncr  10^  o'  . (7) 


TABLE  II.  Surface  tension  for  certain  liquid  metal -organic 
compound  and  liquid  metal-air  interfaces;  also  estimates  of  the 
critical  point  temperatures  for  these  metals  (at  which  the  surface 
tension  goes  to  zero).  The  fourth  column  gives  the  calculated 
critical  reduced  wavelength  = 7 \/?rr  (transition  from  stability 

for  shorter  wavelengths  to  instability  for  longer  wavelengths) 
for  a small  amplitude  disturbance  in  a molten  metal  sheet  10“^cm 
thick  of  laminar  construction  with  N = ICr  lamina,  under  conditions 
where  the  acceleration  is  rj  cm/sec2.  The  column  ^fastest  = 

3a  ^crit  gives  the  reduced  wavelength  for  which  growth 

small  amplitude  theory  predicts  the  greatest  grgwth  constant,  i 

c<.  = o CrmAX  = (maximum  value  of  (gk  - <rA3/>°  )2  = (2g  / 3 Xf.astest)2* 
This  growth  constant  itself  is  tabulated  in  the  next  coluim.^pwfcn 
Consider  the  initial  surface,  prepared  however  so  well,  and 
Fourier  analyze  the  irregularities  in  this  surface,  and  enquire 
as  to  the  typical  amplitude  associated  with  disturbances  of  the 
reduced  wavelength  ^fastest  * Let  i-fc  be  re<luired>  for  the  saka 

growth 

of  illustration,  that  these  disturbances  not  grow  duringatime 


Interval  of  the  order  of  1M-&* o (for  example)  to  the  point 

(amplitude^  3>.3.  ) where  the  non-linear  effects  of  Taylor 

instability  take  hold.  Then  the  initial  amplitude  A must  be 
less  than  ru  expc-^.^.ti  . This  preposterous 

number  is  tabulated  in  the  last  column.  The  impossibility  to 
attain  any  such  perfection  in  the  surface  shows  the  hopelessness 
of  counting  on  surface  tension  to  inhibit  Taylor  instability  in 
the  present  circumstances. 


Interface 

Torit(°K) 

T(°K) 

2 

<r(erg/cm  ) 

^2rit(cm) 

sec-D 

Amax(cm) 

allow- 

■■Hi 

able 

Mercury- 

Benzene 

2 93 

3 57 

1.7  Xto'^ 

-y 

X to 

x,*8 

.J  -£ 

Mercury- 

Air 

> /S‘23 

2 88 

987 

-y 

— 3.3  X /« 

xio8 

-V  -1 

^3.3A*>  X€ 

Copper^ 

IH<>9 

II  o 3 

-v3.  <x/o"^ 

^(,.2  x /o'* 

xl°8 

. J -!• 

x*  Ve 

Aluminum- 

Air 

• 

973 

89o 

-y 

-1 Jlo.z  X/o 

~8-l  xio7 
1 

-3  -S 
x/o  xe 

It  is  evident  from  the  numbers  for  ^c*it  in  the  table  cal- 
culated for  a number  of  laminations  as  great  as  N = 10^  that  it 
has  not  been  of  much  use  to  laminate  the  metal  sheet  to  inhibit 
Taylor  instability!  In  addition,  the  metal  sheet  is  thought  of 
as  driven  by  magnetic  fields  so  great  and  comes  to  temperatures 
so  high  that  it  is  carried  from  the  liquid  state  past  the  critical 
temperature — where  the  surface  tension  vanishes — into  the  gaseous 
state  where  the  normal  concept  of  surface  tension  has  to  be 
abandoned . 

The  natural  scale  of  dimensions  at  the  beginning  of  the  merger 


process  - as  distinguished  from  the  scale  of  fabricational  imper- 
fections ( which  are  overlooked  here,  see  section  IV)  - is  so  small 
(Table  I)  and  the  accelerations  so  high,  10^  (or  lO1^)  cm/sec2 
that  it  can  well  be  asked  (Table  III)  whether  viscosity  may  not 
be  more  important  than  surface  tension  in  governing  the  precise 
value  of  this  natural  scale  of  sizes. 


TABLE  III. 


Critical  acceleration  Ecrit  = 


Here  o'  is  the  surface  tension  (erg/cm2),  J>  is  the  density 
(g/enp)  and  ~y-  VIP  is  the  kinematic  viscosity  (cm2/sec). 

For  accelerations  much  higher  than  gcrit,  viscosity  dominates 

over  surface  tension  in  determining  the  circular  wave  number  At-?. 
and  wavelength  sir/Jtf.g.  of  the  fastest  growing  small  amplitude 

disturbances  and  the  exponential  growth  constant  (sec1 ) 

of  these  Taylor  Instabilities.  For  accelerations  less  than 

the  surface  tension  dominates  over  viscosity  in  determining  the 
circular  wa^e  number  and  the  growth  constant  of  the  small  amplitude 
disturbance.,  which  multiply  at  the  highest  rate: 
and  ^ ( f3P/o'y,‘/  . When  surface  tension  dominates  it 

stabilizes  small  amplitude  disturbances  whose  reduced  wavelengths 
are  of  the  order  of  the  reduced  wavelength  of  the  fastest  grow- 
ing disturbances/  jfcy.y.  , and  less.  When  viscosity  dominates  it 
does  not  stabilize  small  amplitude  disturbances  of  any  wavelength, 
but  it  does  reduce  the  rate  of  multiplication  of  disturbances 
with  3t<  /t/.j..  relative  to  what  the  growth  constant  would  be  in 
the  absence  of  viscosity,  . Each  material  in  the 

table  is  considered  under  circumstances  where.it  is  pushed  by  a 
medium  (gas  or  magnetic  field  or  radiation)  of  density  negligible 
in  comparison  with  the  density  of  the  pushed  medium.  The  temp- 
eratures in  the  table  can  be  converted  to  electron  volts  by 
dividing  the  conversion  factor  1.15  x 10^  °K/electron  volt. 


Pushed  Medium 


T(°K) 


293 
1500 
lpOOD 
2x10/ 
dependent 
upon  cir- 
cumstances 


^CRlT 


The  viscosity  listed  here  for  fully  ionized  hydrogen  under 
conditions  of  interstellar  density  refers  to  the  conditions  of 
ideal  laminar  flow  such  as  are  assumed  in  conventional  definitions 
of  viscosity.  In  other  words  z/3"?//9  has  been  estimated  from  the 
formula (V  ) : 

if  or™*'T),'i/Tr<y-tS>  , T is  the  temperature  (°K), 
k the  Boltzman  constant,  Tr  is  the  collision  cross  section 
for  a rigid  spherical  molecule,  D is  the  coefficient  of  diffusion. 

In  actuality,  the  interstellar  medium,  like  the  earth’s 
atmosphere,  will  typically  be  ir  turbulent  motion  already  before 
it  is  subject  to  pressure  from  magnetic  fields  or  radiation  such 
as  can  generate  Taylor  instability.  Under  these  conditions, 
the  exchange  of  momentum  between  one  part  of  the  gas  and  another 
is  brought  about  very  little  by  kinetic  theory  motion  of  individual 
particles,  such  as  are  considered  in  the  elementary  theory  of 
viscosity.  The  momentum  exchange  comes  about  to  an  enormously 
greater  extent  through  eddies  and  vortices  that  carry  whole  masses 
of  the  medium  from  a layer  with  one  average  velocity  to  a layer 
with  another  average  velocity.  The  effective  viscosity  to  be 
used  under  such. conditions  is  discussed,  for  example,  in  Landau 
and  LifshitzC10/,  and  Hirschfelder,  Curtiss  and  Bird*-?'. 


Viscosity,  like  surface  tension,  decrease  * at  first  as  the 
temperature  of  a molten  metal  is  raised  toward  the  critical 
temperature.  However,  contrary  to  the  surface  tension  which  goes 
to  zero  at  the  critical  temperature,  the  viscosity  goes  only  to 
a reduced  value  somewhere  in  the  neighborhood  of  the  critical 
r.emperature  and  pressure  provided  that  the  path  pursued  in  the 
PT-plane  leads  near  the  critical  point.  As  the  temperature  is 
increased  to  a figure  much  higher  than  the  critical  temperature 
the  viscosity  of  the  metal,  like  the  viscosity  of  ether  gases, 
approaches  Maxwell's  law  of  proportionality  to  the  square  root 
of  the  temperature  T (here  and  in  the  following  expressed 
directly  in  energy  units  to  eliminate  the  purely  conventional 
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Boltzmann  factor  k and  its  possible  confusion  with  the  wave 
number  k): 


7 ^ 2 


Mi  aTi* 


C TVn/»^s  , t • 


(8) 


i &7KANS  , I 

Here  is  the  mass  of  an  ion  (i  = l)  or  electron  (i  = 2), 
vi  is  a velocity  typical  of  thermal  agitation  and  <Tj.rans  is  an 
effective  cross  section  for  transport  such  as  is  defined  in 
diffusion  theory ^).  For  very  high  temperatures  it  is  no  longer 
possible  to  treat  this  transport  cross  section  as  independent  of 
velocity  - and  therefore  as  independent  of  temperature . Its 
fall  with  temperature  makes  the  viscosity  rise  with  a power  of 
T higher  than|-.  Table  IV  gives  a very  rough  estimate  of  the 
viscosity  of  copper  at  a temperature  T ^ 1.3  ev  (150C0  °K). 

TABLE  IV.  Order  of  magnitude  estimate  of  contributions  to 
the  viscosity  of  copper  at  T ^ 1.3  ev.  It  is  assumed  that  the 
typical  atom  at  this  temperature  has  been  dissociated  into  one 
Cu^+ion  and  two  electrons. 


Momentum 

1 

- 

1 

| cr 

transport 

Contribution 

Particle 

JL 

^(M^T)2  (gcm/sec) 

(cm2) 

to  y (g/cm  sec) 

Electron 

-^0.5  x IQ"19 

-17 

^ 10 

— 0.005 

Cu++  ion 

^150  x 10"19 

io-16 

0.15 

Total  y ^ 0.16  g/cm  sec 

Kinematic  viscosity  (at  y0 /~'8g/cm2)  estimated  as  2/  = '///)~'0,02 

cm2/sec 
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Consider  a sheet  of  copper  vapor  with  the  kinematic  viscosity 
2 x 10  cm2/«p  sub ject  to  acceleration  by  a pushing  fluid  of 
negligible  density  at  a rate  of  10^  cm/sec  . Then  in  the 


small  amplitude  theory  of  Taylor  instability  those  Fourier  com 
ponents  are  calculated  to  grow  most  rapidly  which  have  a wave 
number  of  the  order 


x/3 


'£a*T0SJ 


• ^ i 2 

/O  0*^1 

y K 10-*  Co^V /<*</, 


i/3 


r - 1 

3 x to  Con 


(9) 


They  have  a calculated  exponential  growth  constant 


•m^jL 


( ?*/*) 


1/3 


/0  Cs>+-,  t / ■**<■ 

2 X lo  * Cs**1 


a/3 

^ i.i  x iQ** (10) 


Here  as  in  the  case  of  surface  tension,  the  scale  of  wave- 
lengths at  which  the  inhibiting  effects  are  significant  are  far 
too  short  to  make  viscosity  of  any  relevance  at  the  scale  of  in- 
stabilities which  are  of  concern. 


III.  PRIMITIVE  STATISTICAL  ANALYSIS  OF  REGIME  OF  BUBBLE 

AMALGAMATION 

So  far  it  has  been  enough  to  consider  the  rise  of  a 
single  bubble  because  the  starting  assumption  was  complete  periodicity 
with  all  bubbles  identical.  In  actuality  there  will  be  small  differ- 
ences in  size  from  bubble  to  bubble  due  to  slight  departures  from 
ideality  at  the  start  of  the  acceleration  of  the  interface.  These 
differences  will  occur  whether  natural  causes  or  machining  determines 
the  dominant  wavelength.  In  the  one  case,  the  spacing  is  determined 
by  the  circumstance  that  the  growth  constant  has  a smooth  maximum 
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as  a function  of  the  wavelength.  Therefore,  even  a few  extra 
atoms  here  or  there  on  the  interface  will  favor  at  the  start  a 
scale  of  dimensions  which  here  is  a little  less  than  the  size  for 
optimum  growth  and  there  a little  greater.  Similarly,  in  the 
case  where  machinery  or  other  techniques  have  built  into  the 
initial  surface  irregularities  with  a unique  characteristic  length, 
there  will  nevertheless  again  occur  small  variations  above  and 
below  this  unique  size  from  one  incipient  bubble  to  another  due 
to  the  presence  or  absence  of  a few  more  atoms  here  or  there  if 
due  to  no  more  mundane  cause  I 

A small  fractional  difference  in  size 

£ = (R2  - Rx)  /\  (11) 

between  one  growing  bubble  and  another  will  cause  a corresponding 
fractional  difference  S/£  hi  the  rate  of  rise  rji  ^ of  the 

two  bubbles  when  they  are  fully  developed.  In  consequence  of  this 
difference  in  rate  of  rise,  one  bubble  will  get  ahead  of  Its 
neighbor  and  dominate  it.  The  two  bubbles  will  come  to  constitute 
in  effect  a single  bubble  with  some  minor  irregularities  in  shape 
(Pig.  3)*  This  process  of  merger  of  two  bubbles  is  a decisive 
feature  of  the  later  stages  of  Taylor  Instability.  Clear  evidence 
of  repeated  mergers  is  seen  in  typical  photographs  of  Taylor  in- 
stability, such  as  those  taken  by  Lewis  (?)  and  Allred  and  Blount 
It  is  necessary  to  examine  in  more  detail  the  conditions  for 
merger  and  the  consequences  of  merger  in  order  to  arrive  at  any- 
thing approaching  a rational  estimate  of  the  time  required  for 
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Fig . 3 . Factors  governing  bubble  merger  in  idealized  model . Upper 
left:  Representative  bubble  pattern  as  viewed  from  aboce.  The 

cross  sectional  area  associated  with  one  bubble  in  this  diagram  is 
identified  with  TTR2  to  define  the  effective  radius  R.  Upper  right: 
Bubble  or  finger  No.  39  is  idealized  as  engulfing  finger  No.  35  as 
soon  as  the  effective  height  z on  of  the  larger  one  exceeds  the 
effective  height  of  the  smaller^one  by  the  radius  of  the  smaller 
one;  that  is,  as  soon  as  Z39  (point  El  ) comes  to  equality  with 
z35  + R35  (point  0 ) . The  arrows  recall  that  the  rate  of  rise 
of  the  larger  bubble  is  greater  than  the  rate  of  rise  of  the 
smaller  finger,  a§  is  required  for  capture  of  the  smaller  one: 
dz/dt  = 0.48  (gR)a.  Instead  of  the  tyue  time  t it  is  convenient 
to  use  the  "reduced  tiqie",  X = 0.48  g2t,  so  that  the  race  of 
rise  becomes  dz/dt  = Ri.  Lower  diagram:  Graphical  procedure  to 

determine  which  pair  of  bubbles  merge  first.  For  simplicity, 
attention  is  limited  to  the  system  of  eight  fingers  outlined  by 
a heavy  line  in  the  upper  left  hand  diagram.  The  calculated  value 
of  the  reduced  time  interval  A'C  required  for  No.  39  to  capture 
No.  35  is  determined  by  the  slope  (with  sign  reversed)  of  the 
designated  line  in  the  diagram.  The  line  from  No.  39  to  No.  42 
is  irrelevant  because  it  corresponds  to  two  bubbles  which  are  not 
in  contact.  For  fingers  in  contact,  the  line  with  the  lowest 
slope  corresponds  to  the  earliest  merger.  A new  diagram  of  this 
type  was  constructed  for  each  successive  stage  in  the  process  of 
bubble  merger  in  the  process  of  working  out  the  illustrative 
numbers  in  Table  V. 
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bubbles  to  break  through  a slab  of  liquid  of  finite  thickness. 

No  detailed  hydrodynamic  analysis  of  the  merger  process  is 
available.  Therefore,  order  of  magnitude  considerations  are 
demanded.  For  this  purpose,  the  actual  situation  will  be  replaced 
by  an  idealized  model,  here  called  the  "Maine  bubble  model" 
because  the  first  illustrative  calculations  on  the  model  were 
made  in  Maine.  In  this  model,  each  bubble  will  be  considered  to 
stake  out  for  itself  a certain  area  when  it  is  projected  on  the 
ideal  initial  interface.  This  area  will  be  assumed  to  remain  con- 
stant during  the  time  of  rise  of  the  bubble.  The  area  will  typically 
be  irregular  in  outline.  However,  the  rate  of  rise  will  be  cal- 
culated as  if  the  shape  were  that  of  an  ideal  circle  ;.xth  the  same 
cross  sectional  area  IT  (definition  of  R!).  For  the  rate  of 
rise,  Taylor's  formula  will  be  adopted, 

rs  = (0.48)  (gR)^  (12) 

appropriate  for  the  case  of  a pushing  fluid  of  negligible  density. 
Suitable  changes  can  be  made  to  correct  the  analysis  when  the  push- 
ing fluid  has  a density  appreciable  in  comparison  with  that  of 
the  pushed  fluid. 

The  volume  of  gas  which  fills  the  bubble  above  the  ideal 
initial  surface,  divided  by  the  cross  sectional  area  of  the  bubble, 
defines  an  effective  height  z for  the  bubble  which  is  somewhat 
lower  than  the  vertex. 

Two  bubbles  will  be  assumed  in  the  model  to  merge  into  one 
when  the  larger  one  is  sufficiently  far  ahead  of  the  smaller  one. 
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The  lead  required  for  merger  will  be  very  little  when  the  smaller 
bubble  is  very  small  compared  with  its  winning  neighbor.  The 
necessary  lead  will  be  greater  when  the  effective  radius  of  the 
smaller  bubble  constitutes  a larger  fraction  of  the  radius  of  the 
larger  one.  This  qualitative  aspect  of  the  merger  process  will  be 
sharpened  up  iri  the  model  to  the  following  rule: 

1st  rule.  Two  bubbles  will  merge  when  the  larger  one  of 
them  has  a lead  on  the  smaller  as  great  as  or  greater  than 
the  radius  of  the  smaller  one: 

z^  - z > R-  . (13) 

T — 

The  other  two  rules  adopted  in  the  model  for  the  merger 
process  are: 

2nd  rule.  Conservation  of  cross  sectional  area, 

irft^  = Trfl*  +irni  . (14) 

This  rule  provides  a way  to  calculate  the  effective  radius, 

R of  the  merged  bubble, 
m 

3rd  rule.  Conservation  of  volume, 

TT  R2m  zm  « TTr|  z+  +7TR?  z_  . (15) 

From  rules  2 and  3 follows  an  expression  for  the  effective 
height  of  the  new  bubble  at  th^noment  of  formation: 

z = R+  z+  + R-  z- 

m — • (16) 

As  starting  cona.  tions  for  the  model  one  has  to  specify  the 
quantities  R and  z for  each  bubble  at  the  time  t = o and  has  to 
give,  in  addition,  the  connectivity  of  the  bubble  pattern,  either 
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by  way  of  a qualitative  diagram  like  Pig.  3>  or  by  way  of  a listing 
of  all  bubble -to -bubble  contacts  which  exist  at  the  time  t = o. 

In  the  example  in  Pig.  3,  these  connections  are  1-2,  1-22,  1-21, 
1-20;  2-3,  2-23,  2-22;  « * . ; 37-38. 

Table  V illustrates  the  evolution  in  time  of  an  arbitrary 
initial  pattern  of  eight  bubbles  calculated  according  to  the  rules 
of  the  foregoing  model. 


TABLE  V.  Development  of  a sample  pattern  of  eight  bubbles 
(Fig.  3)  as  computed  from  rules  1,2,3  of  the  Maine  bubble  model. 

The  initial  conditions  were  arbitrarily  adopted  to  correspond  to 
(1  ) a fractional  variation  of  the  order  of  50  percent  in  bubble 
cross  sections  and  (2)  an  Initial  elevation  for  each  bubble  equal 
to  the  effective  radius  of  that  bubble.  This  condition  is  meant 
to  correspond  qualitatively  to  the  circumstance  that  the  larger 
disturbances  have  typically  risen  to  a greater  height  than  have 
the  smaller  ones  at  the  stage  when  the  bubbles  have  come  into 
the  non-linear  regime  of  linear  rise  with  time.  However,  the  here 
assumed  proportionality  of  2L  to  the  first  power  of  R.  has  no  basis 
in  theory  nor  any  special  significance.  The  assumption  is  made 
only  to  give  definite  and  simple  starting  conditions  for  this 
particular  problem.  In  the  table,  the  quantity  "reduced  time" 
x is  an  abbreviation  for  the  quantity  (2/3)  g^t.  In  tegms 
of  x the  formula  for  the  velocity  of  rise  becomes  dz/dr  = R'2\ 

At  each  stage  of  the  calculation,  every  contacting  pair  of  bubbles 
is  considered  in  turn,  and  for  each  such  pair  the  reduced  time 
to  merger  is  calculated  in  each  phase  of  development  that  one 
merger  is  listed  which  is  calculated  to  occur  earliest. 


Bubble 

26 

33 

34 

35 

39  40 

4l 

42 

Reduced  time 

R 

0.3l 

1.21 

1.10 

0.74 

1.44  1.09 

0.64 

0.72 

for  merger 

'Z 

0.81 

1.21 

1.10 

0.74 

1.44  1.09 

0.64 

0.72 

AX,  X = 2 AX 

dz/dr 

0.90 

1.10 

1.05 

0.86 

1.20  1.04 

0.80 

0.85 

Bubble 

26 

33 

34 

35-39 

Up. 

Ml 

42 

R 

0.8l 

1.21 

1.10 

1.62 

1.09 

0.64 

0.72 

Z 

0.8l 

1.21 

1.10 

1.29 

1.09 

0.64 

0.72 

= 0.23 

dz/dr 

0.90 

1.10 

1.05 

1.27 

1.04 

0.80 

O.85 

X = 

0.23 
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Bubble 

26  33-41  34 

35-39 

40 

42 

R 

0.81  1.37  1.10 

1.62 

1.09 

6.72 

Z 

0.8l  1.08  1.10 

1.29 

1.09 

0.72 

0.89 

dz/dr 

0.90  1.17  1.05 

1.27 

1.04 

0.85 

rr  = 1.12 

Bubble 

26-35-39  33-41  34 

40 

42 

R 

1 . 8l  1.37  1.10 

1.09 

0.72 

Z 

0.80  1.08  1. 

10 

I.09 

0.72 

1.13 

dz/dr 

1.34  1.17  1.' 

05 

1.04 

0.85 

'C  =2.25 

Bubble 

26-35-39  33-41-42 

34 

40 

R 

l.8l  1.55 

0.80  1.00 

1.10 

1.09 

AK~~=  1.27 

Z 

dz/dr 

1.10 

1.09 

r = 3.52 

1.34  1.24 

1.05 

1.04 

Bubble 

26-35-39-40  33-41 

-42 

34 

R 

2.12  1,55 

1.10 

Ar^=  3.21 

Z 

dz/dr 

0.88  1,00 

1.10 

**  = 6.73 

1.46  1.24 

1.05 

Bubble 

26-34-35-39-40 

33-41- 

42 

R 

2.38 

1.55 

ZVC^=  5.40 

Z 

0.93 

1.00 

r = 12.13 

dz/dr 

1.54 

1.24 

Bubble 

26-34-35-39-40-33- 

cvj 

■=r 

1 

H 

R 

2.84 

Z 

0.95 

dz/dr 

1.68 
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The  average  velocity  of  the  system  of  fingers  increases  step 
by  step  as  the  mergers  proceed  until  finally  the  entire  cross 
sectional  area  is  occupied  by  a single  finger  rising  at  a constant 
velocity.  Here  it  is  assumed  that  the  slab  of  fluid  is  thick 
enough  so  that  breakthrough  does  not  occur  until  finger  amalgama- 
tion is  complete.  A thinner  slab  of  heavy  fluid  will  be  penetrated 
by  a finger  of  the  lighter  fluid  which  is  still  not  merged  with  all 
other  fingers. 

It  follows  from  this  analysis  that  a slab  of  fluid  very  thick 
irykomparison  with  its  effective  radius  will  almost  always  be 
penetrated  by  a single  finger.  Independently  of  the  original  pattern 
of  the  fingers.  Moreover,  the  time  for  merger  of  the  fingers  in 
this  case  will  constitute  only  a small  fraction  of  the  much 
longer  time  taken  for  the  final  big  rapidly  moving  finger  to  break 
through.  Therefore,  an  estimate  of  the  time  for  breakthrough  can 
be  made  in  this  case,  from  the  formula 


/ 


breakthrough  time 
for  a slab  of  fluid 
thick  in  comparison 
with  its  effective 
radius 


(thickness  of  slab) 
Ptf; 


(2/ 

(17) 
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In  this  time  the  slab  has  advanced  through  the  distance 


/ distance  \ ^ 
\ of  advance) 


(l/2)  (acceleration)  (time  )£ 


(9/8) (thickness  of  slab)  /(effective  radius)^ 

(18) 
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a quantity  which  is  large  when  compared  either  with  the  effective 
radius  or  with  the  thickness  of  the  slab. 

An  estimate  of  the  breakthrough  time  in  the  case  of  a thinner 
slab  demands  an  analysis  of  the  details  of  merger.  In  the  case 
of  the  eight  fingers  considered  in  Table  V,  the  velocities  change 
with  time  as  indicated  in  Fig.  4. 

The  spread  in  velocities  throughout  most  of  the  stage  of 
cannibalization  is  roughly  of  the  order  of  one  quarter  to  one 
half  of  the  average  velocity  itself  at  each  instant.  This  average 
velocity  increases  in  a stepwise  fashion  with  time.  However,  if 
the  increase  is  represented  in  idealized  form  by  a curve,  then  a 
straight  line  curve  is  not  out  of  place  to  describe  the  results 
of  the  computations.  The  rate  of  increase  of  velocity  with  time 
read  off  of  such  a curve  lies  between  0.02  g and  0.05  g in  the 
example . 

Abstracting  from  this  special  example  one  is  led  to  consider 
a pattern  of  fingers  which  develops  in  time  according  to  a similarity 
on  scaling  transformation,  as  follows: 

(l)  An  average  rate  of  rise  for  the  fingers  which  are 
present  at  any  given  time  which  is  given  by  the  expression 
v = (average  velocity)  = k^  gt, 

where  k is  a small  dimensionless  numerical  factor  of  the 
1 

order  of 

0.02  to  0.05.  (19) 
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Pig.  4,  Velocity  as  a function  of  time  for  the  fingers  listed  in 
Table  V.  Each  horizontal  line  corresponds  to  one  finger.  The 
line  starts  when  the  finger  is  formed  and  disappears  when  the  finger 
merges  with  another.  Some  foreshadowing  can  be  seen  in  this  diagram 
of  the  central  points  of  a statistical  analysis: 

(1 ) A statistical  distribution  in  rise  velocities  at  any 
given  time.  This  distribution  in  the  example,  neither  clearly 
narrows  up  as  time  advances  nor  clearly  broadens  out  when  expressed 
in  terms  of  the  fractional  variation  of  the  velocity  about  the  av- 
erage velocity  which  obtains  at  that  instant. 

(2)  A roughly  linear  increase  of  this  average  rise  velocity 
with  time,  leading  to  the  concept  of  an  effective  acceleration 
associated  with  bubble  rise  during  the  regime  of  continued  merger. 

The  two  dashed  lines  correspond  to  values  of  this  effective  accelera- 
tion of  0.02  g and  0.05  which,  therefore,  serve  in  the  particular 
example- as  lower  and  upper  limits  for  this  effective  acceleration. 


(20) 


(2)  An  average  effective  finger  radius 

^ I oR-  Kav=A:^ 

v;here  k is  another  dimensionless  factor. 

2 

(3)  A statistical  distribution  of  finger  radii  about  the 
average  radius  given  by  something  qualitatively  like  a 
Gaussian  curve,  with  a fractional  spread  in  radii  which  does 
not  change  in  time  as  the.  fingers  merge: 

[ (R-  Rav)2]av  /F£v  = <f 2 = COHSTANT  . 

If  a distribution  as  special  as  the  Gaussian  distribution 

were  to  be  assumed,  one  would  have  the  formula 

/ number  of  fingers  ^ / total  number  of^  dR  f (R-Rav  )2 

d f with  radii  between)  = fingers  at  the  / (2  7T)i  R^  6xp /"  0 tn  r ^ 
y R and  R+dR  j \ given  time  J l d b J 

(21) 

not  only  for  the  radii  themselves,  but  also  for  the  cross 
section  area.  For  the  rise  velocity  and  for  the  effective 
height  distributions  are  assumed — expressed  in  fractional 
deviations  from  the  average — which  do  not  change  with  time. 

This  statistical  assumption  is  taken  to  apply  throughout 
the  regime  of  continuing  mergers,  apart  from  the  first  few 
generations  of  mergers  where  the  specialities  of  the  start- 
ing mechanism  show  themselves,  an^part  from  the  final 
stage  (if  breakthrough  has  not  already  taken  place)  where 
only  a few  fingers  are  left. 

It  is  possible  to  make  these  assumptions  plausible  on  dimen- 
sional grounds.  In  the  regime  under  consideration,  the  fingers 
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have  grown  in  size  to  dimensions  so  great  compared  to  the  characteris- 
tic 3:  set  by  surface  tension  and  acceleration  (or  by  viscosity 

and  acceleration)  that  surface  tension  (and  viscosity)  have 
nothing  to  do  with  the  main  features  of  the  rise  phenomenon.  At 
the  same  time  the  fingers  are  small  enough  in  comparison  with  the 
dimensions  of  the  slab  of  fluid  that  these  macroscopic  dimensions 
also  can  have  no  significant  effect  on  what  is  going  on.  Therefore, 
the  only  physical  quantities  left  to  dominate  the  situation  are 
the  acceleration  g and  the  elapsed  time  t.  Out  of  these  magnitudes 
there  is  only  one  way  to  construct  a quantity  with  the  dimensions 
of  an  average  finger  radius  (Eq.  19)  and  a quantity  with  the  dimensions 
of  an  average  rise  velocity  (Eq.  20). 

These  dimensional  arguments  do  not  justify  formulas  of  the 

type  (19)  and  (20 ) with  universal  coefficients  k and  k , indepen- 

2 

dent  of  the  original  spread  in  finger  sizes.  It  is  even  possible 
to  point  to  a perfectly  conceivable  situation  in  which  the  dis- 
tribution in  finger  sizes  deviates  enormously  from  a Gaussian 
curve,  not  only  in  the  first  few  generations  of  mergers,  but  also 
right  up  to  the  stage  where  only  a few  bubbles  are  left.  (Pig.  5). 
Nevertheless,  even  in  this  extreme  case  a simple  analysis  makes  it 
reasonable  to  think  of  the  effective  rate  of  rise  as  being  unaffected 
by  the  anomaly  in  the  distribution  of  finger  sizes. 

Rather  than  the  big  finger  in  Pig.  5 helping  the  little  ones 
to  move  forward  faster,  the  little  fingers  hold  back  the  big  one, 
if  one  can  take  the  Maine  model  as  a guide  to  the  true  state  of 


LIGHT  FLUID 


Pig.  5.  Case  of  anomalous  distribution  in  finger  sizes.  The 
larger  finger  rises  much  more  rapidly  than  the  small  ones  and 
engulfs  them  one  by  one  until  it  spans  the  entire  cross  sectional 
area  of  the  slab  of  fluid. 
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affairs.  According  to  one  of  the  rules  of  that  model,  a large 

bubble  of  effective  height  z engulfs  smaller  ones  of  average 

effective  height  z as  soon  as  the  difference  z - zav  rises  to  a 

value  as  great  as  the  average  effective  radius  R of  these 

av 

fingers.  The  consequence  is  this,  that  the  big  finger  grows  in 
radius  faster  than  the  little  ones,  but  does  not  get  out  of  step 
with  them  in  vertical  elevation.  In  each  interval  of  time  dt  it 
eats  whatever  number  of  little  bubbles, dN,  is  required  to  dilute 
its  rise  back  to  parity  with  its  fellows.  The  new  effective  height 
is  found  by  adding  the  volume  of  the  old  finger  and  the  newly 
assimilated  ones  and  dividing  by  the  cross-sectional  area  of  this 
collection.  If  too  many  small  bubbles  have  been  taken  in  during 
the  interval  dt,  this  effective  height  will  be  reduced  too  much 
and  new  cannibalization  cannot  go  on  until  the  faster  rate  of  rise 
of  the  larger  bubble  once  more  allows  it  to  get  ahead.  The  converse 
also  applies — if  the  big  bubble  assimilates  too  few  small  ones,  its 
elevation  will  not  be  diluted  back  enough  to  make  up  for  its  faster 
rate  of  rise.  Therefore,  it  will  get  ahead  and  in  the  next  time 
interval  take  in  a more  than  normal  number  of  fingerlets.  Thus, 
in  the  model,  there  is  a stability  about  the  mechanism  through  which 
the  little  bubbles  govern  the  rate  of  rise  of  the  big  one. 

This  reasoning  can  be  expressed  in  mathematical  terms.  Let 
the  average  rate  of  rise  of  the  small  fingers  be  expressed  in  the 
form  (19)  expected  from  similarity  arguments, 

v (t)  = lc j gt  (22) 

av  * 
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and  let  their  average  height  and  average  radius  be  written  in  a 
form  consistent  with  this  expression: 


V = (*j.  A)  St2 
r (t)  - (^  /o.m2  gt2 


(23) 

(2ft) 


The  big  finger  is  always  at  the  point  where  the  smallest  addition 
to  its  height  makes the  difference  between  more  cannibalization  and 


none  at  all;  or,  according  to  the  model. 


* 

z - z 


av  - Rav 


Thus  one  finds 


z~  (t)  = [(kj/2)  + (ltj/0.48)2]  gt2 


(25) 


(26) 


The  rate  of  rise  of  the  big  finger  calculated  from  this  expression 
is 

v*(t)  = C.dz  */dt  = [k  -i-  (k2/0.12)]  gt  # (27) 

Any  value  of  the  effective  radius  R of  the  large  bubble  is  com- 

* 

patible  with  this  expression.  However,  the  larger  R happens 
to  be,  the  greater  is  the  natural  rate  of  rise  associated  with  a 
finger  so  big, 

/ natural  rate  of  \ « i 

\ rise  of  big  finger  / = (0.48)(gR  )?  (28) 

and  the  more  little  bubbles  it  must  assimilate  per  second  to  hold 

its  rate  of  rise  back  to  the  figure  (27).  The  "drop  back"  due 

* 

to  eating  out  from  the  effective  radius  R to  the  effective  radius 


* * 

R + dR  in  the  time  dt  is 


,*r.*2 


* . 


2 /!«/£  -*-4ifc»  2 -vTC^tAY  / '*  '• 5 ;C  - Z - 


‘ * z R"c  + zav  2R*  dR* 


*p  # * 

R + 2R  dR 


= (2dR*/R*)(z*-  z ) 


= (2dR*/R*)  R. 


av' 


av 


(29) 


The  rate  of  drop  back  on  account  of  assimilating  smaller  bubbles 
has  to  satisfy  in  the  model,  the  equation 


rate  of  drop  back  \ /natural  rate  of 

due  to  assimilation  J = rise  of  big  finger  ) - I required  to  remain  ir 
of  smaller  bubbles  J v in  absence  of  growth/  ^ step  with  smaller 

bubbles  > 


\ f slower  rate  of  rise 


(30) 


or 


(2  Rav/R  )(dR  /dt)  = (0.48)(gR>  - (0.48) (gRav)2  . (31) 

Under  the  conditions  contemplated  in  Fig.  5,  the  radius  R*  is 
large  compared  to  Ray.  Therefore,  it  is  appropriate  to  neglect 
the  last  term  in  (31).  Then  Eq.  (31)  for  the  growth  of  the  big 
finger  may  be  integrated  to  give  the  effective  radius 

•>2. 

(32) 


R (t)  = R*(to) 


Ik, 


r 


l ?Jt,  4 [ R*(t*)/o.S£  R^W]-[R  (tv /o-sre  j 


Here  the  radius  R (t)  of  the  small  fingers  has  been  taken  to 

civ 


follow  Eq.  (24). 

To  fix  ideas,  assume,  for  example 
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( advance  of  bubbles  } 
l relative  to  liquid/  ^ " 

.Liquid  / 

13  . , 2 


k = /advance  of  \ 

1 [slab  of  liquid  / 


breakthrough  \ 

constant"  J = 0.02; 


g = 2 x 10 


cm/sec  ; 


(dimension  fixed  \ 

by  small  amplitude  \ i 

theory  allowing  j ^ 3*8317  (3  /g Pi 

for  surface  tension/ 


/ 


^ 6.6 


(800g/sec  ) 


•I  2 


-5 


= 1.3  x 10  cm; 


L (2  x 10^  cm/sec2) (10  g/cm3)  J 

vav(tQ)  - (0.48)(gRav)i  = 7-7  X 10  3 


cm; 


zav(tQ)  = (0.23/2kx)  Rav  (Eqs.  23,  24)  = 7.8  x 10"5 


cm. 


(33) 


These  "initial"  values  correspond  to  a time  coordinate 

t = 2.0  x 10“®  sec  in  Eqs.  (23)  and  (24).  At  any  later  t the 

calculated  effective  radius  of  the  small  fingers  is 

10  9 P 

Rav(t)  = (3-3  x 10  cm/sec  ) t . (34) 


Assume  that  a flaw  somewhere  in  Lhe  original  surface  gives  rise 
to  a single  disturbance — in  the  midst  of  these  smaller  bubbles — 


with  an  effective  radius 

R (tQ)  = 10-3  cm. 


(35) 


Then  the  calculated  radius  jf  this  dominating  finger  at  any  later 
time  is 


R*(t  ),  = 10"3 


cm 


0.18 

0.18  + 2 x 10“ 3 cm  - 2 x 10~^cm 


(36) 


Rav  (^ ) 


1.3  x 10”5cm' 


v. 


The  denominator  vanishes  and  the  finger  radius  goes  to  infinity 
when  Rav  has  increased  by  the  very  small  amount 

ARav  = (0.18)2  (1.3  x 10-5  cm)2  /4  x 10_3cm  = 1.4  x l(T9cm 
to  the  value 

Rav(nev;)  = Rav  (tQ)  + /JR  = 1.3x  10~5  + 1.4  x 10"9cm.  (37) 

At  this  instant,  according  to  the  model,  the  big  finger  has  grown 
so  much  that  it  engulfs  the  entire  fluid.  Until  this  instant, 
the  little  ones  govern  the  rate  of  rise;  thereafter,  the  big  one 

does,  and  the  velocity  keeps  the  constant  value 

. . , > - < 

sRftyXviw*/*  *')>-  • 

trike  J 

=(0.48)  g’2"(area  of  slab/TT'  )4.  (38) 

The  example  just  discussed  would  make  it  appear  that  one 

finger  substantially  larger  than  its  fellow^jill  quickly  come  to 

dominate  the  flow  pattern  of  the  entire  mass  of  fluid.  How  much 

validity  should  be  ascribed  to  this  consequence  of  a rather  idealized 

model  of  bubble  merger?  (l).  Available  experimental  information 

(Davies  and  Taylor^;  Lewis Allred  and  Blount  Chang^11)  ) 

does  not  show  any  such  sudden  rise  of  one  finger  to  dominance. 

However,  these  experiments  are  not  relevant  anyway  to  the  point  at 

issue.  The  original  size  of  all  fingery-and  the  thickness  of  the 

slab  of  fluid — were  not  great  enough  compared  to  the  characteris- 

x 

tic  dimension  6.6(  <r/  gj>)2  set  by  surface  tension  and  by  the 
acceleration.  Therefore,  surface  tension  greatly  changed  the 


33 


phenomena  which  inertia  and  pressure  would  otherwise  have  brought 
into  evidence — and  which  are  the  object  of  concern  here. 

(2)  A qualitative  look  at  the  hydrodynamics  of  bubble  merger 
in  the  absence  of  surface  tension  would  suggest  that  a finite  time 
is  required  to  complete  the  process  of  bubble  merger.  The  debris 
from  the  smaller  finger  has  to  be  washed  down  the  surface  of  the 
newly  enlarged  finger  before  this  object  can  be  regarded  as  a unit. 
Only  then  can  it  be  expected  to  live  up  to  its  new  cross-sectional 
dimensions  in  its  ability  to  engulf  a new  bubble.  Consequently, 
one  would  expect  a certain  limitation  on  the  rate  at  which  a 
larger  bubble  can  spread  its  Influence  sideways  to  take  over  one 
smaller  bubble  after  another.  No  such  limitation  of  rate  was 
built  into  the  Maine  model  of  finger  merger.  Therefore,  it  is 
permissible  to  doubt  the  prediction  of  the  model  that  a significant- 
ly larger  thatj&verage  bubble  will  take  over  the  whole  fluid 

before  the  other  fingers  have  had  time  even  to  double  in  height. 

(3)  It  is  conceivable  that  the  distribution  of  finger  sizes 
ii^ayleigh -Taylor  instability  has  a certain  analogy  to  the  dis- 
tribution of  electron  velocities  in  a gas  subject  to  an  electric 
field  in  this  respect, that  there  is  a sharp  distinction  between  a 
statistically  stable  situation  and  a runaway  situation.  In  the 
electrical  case,  the  critical  parameter  is  the  applied  electric 
field.  In  the  hydrodynamic  case,  the  difference  between  (l) 
runaway  and  (2)  approach  to  a standard  distribution  is  governed — 
on  this  view--not  by  the  acceleration  (regarded  as  fixed)  but  by 
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the  initial  distribution  of  finger  sizes  itself.  In  other  words, 
a distribution  which  departs  not  too  greatly  from  the  canonical 
distribution  is  conceived  to  approach  in  time  the  canonical  one. 

The  characteristic  scale  R „ of  this  canonical  distribution  is,  of 
course,  itself  increasing  at  the  same  time.  On  the  other  hand,  a 
distribution  of  finger  sizes  which  departs  too  greatly  (Pig.  5) 
from  the  canonical  one^conceived  to  depart  more  and  more  from 
that  distribution  and  to  lead  to  one  or  more  runaway  bubbles. 

In  summary,  present  information  does  not  suffice  to  distinguish 
between  the  following  three  possibilities: 

(1)  One  bubble  always  quickly  runs  away . 

(2)  Depending  upon  the  degree  of  anomaly  in  the  original 
distribution  of  sizes,  either  one  bubble  runs  away,  or  a 
standard  statistical  distribution  in  sizes  is  attained,  with 
a corresponding  standard  law  of  growth. 

(3)  A standard  distribution  and  rate  of  rise  are  almost  always 
reached . 

Possibilities  with  a light  fluid  to  push  a thin  slab  of  heavy 
fluid  a great  distance  would  seem  to  depend  very  heavily  upon  the 
correctness  of  alternative  (2)  or  (3). 

How  rapid  will  be  the  rate  of  breakthrough  if  (2)  or  (3) 
apply?  The  following  order  of  magnitude  analysis  of  this  question 
is  assumed  for  simplicity: 

(l)  something  like  the  Maine  model  for  bubble  merger  plus 


35 


..  >A.i  .1  V^'W^1  ' ^--*JL-’ j - l.l W ^ LWfcq?j^ 1 f,liJ.^-gr!  %r~  ^.■••^,IW 


(2)  stability  of  a canonical  distribution  of  finger  sizes. 


In  addition,  for  simplicity  and  definiteness,  this  distribu- 
tion will  be  taken  to  follow 

(3)  a Gaussian  distribution  lav;,  with  a parameter  8 of 
fractional  spread  in  radii  equal  to  0.5: 


d (number  of  fingers) 
(total  number  of  fingers) 


dR, 

WtHctt 


exp 


(R-«av) 

2(R  ~S) 


8 = 0.5. 


(39) 


Expression  (39)  predicts  a non-zero  probability  for  fingers  of 
negative  radius.  However,  the  absolute  value  given  for  this 
probability  is  so  lov;  that  this  minor  drawback  of  the  Gaussian 
formula  can  be  disregarded  in  favor  of  mathematical  simplicity. 

Now  to  estimate  the  rate  of  evolution  of  such  a canonical  distribu- 
tlon  assuming  its  self-perpetuating  character! 

Let  R^  ^ Rav  (1-  8/2.)  represent  the  effective  radius  of 
one  bubble  which  is  soon  to  be  eaten.  The  larger  one  which  eats 
it  has  a radius  of  the  order  of  Rn  ^ R (l  + /2).  The 

velocities  of  rise  of  the  two  fingers, 

^0.48  (gRavP'(l-  S/2f  and~0.48  (gR^2  (l  + 8/ 2)?, 


will  differ  by  an  amount 
Av  ^ 


o.48  (gR  ) ( S/e). 

civ 


(4o) 


One  must  get  ahead  of  the  other  by  an  amount  42  of  the  order  Rav 
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to  merge  with  it.  The  time  required  for  this  lead  is  of  the  order 


At  ^ (2/0.48  & 


)(Bav/s)5 


(41) 


In  this  time,  two  bubbles  have  amalgamated  to  give  a single  one 
of  roughly  tv.’ice  the  area.  The  same  effect  has  been  happening 

p 

everywhere  else  through  the  fluid.  Thus,  ^Rav  has  been  in- 
creased by  a factor  2;  and  R&v  has  been  magnified  by  2*;  and  the 
average  velocity  of  rise  has  gone  up  by  a factor  24,  from 
A x i 

0.48  (gR  )2  to  2U  x 0.48  (gR  )2.  (42) 

av7  av 


The  increase  in  vfelocity,  divided  by  the  time  A t required  to 
bring  it  about,  defines  the  effective  acceleration  (see  Pig.  4,  6) 
associated — by  the  statistical  model — with  finger  penetration: 


r>J 


( average  acceleration  \ 
of  fingers  with 
respect  toothe  base  line 
which  would  be  associated 
with  acceleration  of 
\ ideally  flat  interface  y 


nJ  A v/A  t 


(24  -1)  0.48  (gRav)2 
(2/0.48  S )(Rav  / g)* 


= (0.022  <f  ) g 


or 

Sav  = k g ^ O.Olg  . (43) 

This  estimate  of  a breakthrough  coefficient  k^,  even  if  the  ideas 
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Fig.  6.  Increase  of  velocity  with  time  as  a consequence  of  bubble 
Jcco^dlng  Jo  the  Primitive  idealized  statistical  model  of 
thl  5?  ??ch  ?teP  the  velocity  increases iby  a factor  and 

2?®  L engUl °i  the  step  increases  by  a factor  24  compared  to  the 
time  occupied  in  the  previous  step.  Compare  this  regular  in- 
crease of  the  velocity  with  that  shown  in  Fig.  4. 
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behind  it  have  qualitative  validity,  is  evidently  uncertain  by  one 
or  more  powers  of  two. 

The  present  model  of  breakthrough  predicts  that  a slab  of 
heavier  liquid  cannot  be  accelerated  by  gas  to  a distance  greater 
than  (l/k1)/V  (25  or)  50  or  100  or  200  (or  400)  times  its  own 
thickness  without  suffering  finger  penetration.  It  can  be  called 
an  optimistic  model  of  Rayleigh  - Taylor  instability  in  this  sense, 
that  it  is  hard — short  of  drastic  measures,  like  freezing  the 
liquid,  etc. — to  think  of  a way  for  a light  fluid  to  accelerate 
a heavy  fluid  to  a greater  distance.  On  the  other  hand,  rapid 
growth  of  larger  than  normal  fingers  could  well  drastically  reduce 
the  distance  over  which  the  liquid  holds  together.  Optimistic  as 
it  is,  the  model  makes  it  hard  to  see  how  a magnetic  field  could 
ever  propel  a 0.01  cm  slab  of  molten  metal  or  plasma  for  a dis- 
tance as  great  as  10  cm. 

Statistical  mechanics  as  applied  to  gases  teaches  that  the 
number  of  molecules  or  the  energy  of  one  extended  region  is  the 
closer  percentagewise  to  the  energy  or  number  of  molecules  in  an- 
other extended  region  of  the  same  volume,  the  larger  this  volume 
happens  to  be  chosen.  Therefore  on  a surface  covered  over  with 
an  enormous  number  of  small  scale  bubbles  or  fingers -one  might 
be  inclined  to  reason-large  scale  differences  between  one  extended 
region  and  another  can  never  arise.  Therefore  the  statistical 
picture  of  bubble  merger  might  be  considered  to  be  complete 
nonsense.  However,  this  reasoning  would  appear  tacitly  to  assume 


that  there  is  some  conserved  quantity  - like  total  energy,  or 
total  number  of  particles-  which  cannot  be  supplied  in  surplus  to 
one  region  except  at  the  expense  of  another  region.  There  is  no 
such  conserved  quantity  for  bubble  advance.  The  sooner  bubbles 
merge,  the  faster  they  can  go  forward  - and  this  not  at  all  by 
depriving  other  regions  of  the  chance  to  get  ahead.  Bubble 
merger  like  hurricane  growth  is  a typical  example  of  a divergent 
process  in  the  sense  of  Langmuir^12 ^ . 


IV.  HOW  FAR  DO  IMPROVEMENTS  IN  THE  SURFACE  SMOOTHNESS  PAY 
OFF  IN  REDUCED  TAYLOR  INSTABILITY 

The  concept  of  a breakthrough  coefficient  — i£  such 

a concept  is  valid — v/ould  seem  to  provide  some  rational  grounds 

for  deciding  at  what  degree  of  perfection  to  stop  trying  to  improve 

the  original  surface. 

The  disturbances  which  have  evolved  by  merger  out  of  the  in- 
evitable initial  small  scale  perturbations  ( X ^ \ fastest  ) 

growing 

are  taken  in  the  statistical  model  to  follow  the  law  of  growth, 

2 


z = (k  /2)  gt‘ 
av  1 


vav  = ^ gt  = 0.48  (gRav) 
Rav  = (k^/0.48)2  gt2 


(44) 


with  a continually  increasing  scale  of  sizes.  In  the  meantime 
the  irregularities,  if  any,  due  to  machining— here  assumed  to  have 
a scale  *maoh  much  greater  than  * fasteBt  growing  -have 
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themselves  been  developing.  Their  normal  rate  of  increase  will  be 
little  affected  by  the  presence  of  the  merging  and  remerging  smaller 
scale  bubble  pattern  until  the  scale  of  that  pattern,  R&v,  has 
reached  the  scale. 


R = 3.83  A . 
mach  mach 


(^5) 


(of  Section  III)  of  the  machining  irregularities..  The  two  scales 
come  into  concordance,  according  to  the  model,  at  a time  (Eq.  44) 

tcavpWc'i'-D'io  = (0.48/^) (3. 83  *mach/g)*-  (46) 


At  this  time  the  calculated  advance  of  the  fingers  which  have 

resulted  from  merger  is 

z = ( (0.48)2  /2k,  ) R 
av  ' 1 av 

= ( (0.48)2/2ki)  Vch 

= 3.83  ( (0.48)2 /a^) 

( ^22  ^Tmach  or  3.5  A ^ch  ^or  ^1  = 


Three  possibilities  present  themselves  at  this  time: 

(l)  The  disturbances  which  originated  from  machining 
Irregularities  have  already  grown  to  an  amplitude  large 
compared  to  (47 ) . In  this  case  they  have  developed  well 
beyond  the  regime  of  th^mall  amplitude  theory,  as  may  be 
seen  from  the  circumstance  that  in  (47)  zav  is  large  compared 
to  Rav.  Therefore,  it  will  be  expected  that  fingers  have 
developed  out  of  the  Initial  fabricational  imperfections 
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and  even  that  these  fingers  have  made  substantial  progress 

towards  merger  with  one  another.  Evidently  the  instabilities 

from  this  source  are  well  ahead  of  those  that  came  from 

^ . They  dominate,  not  only  at 

A fastest  growing 

t concordance  but  thereafter.  Consequently,  it 

can  be  said  that  the  machining  irregularities  are  decisive 
in  determining  the  time  of  breakthrough  in  case  (l). 

(2)  The  disturbances  which  originate  from  imperfections 

in  the  surface  of  reduced  wave  length  ^ mach  have  Srown 

by  the  time  t of  Eq.  (46)  to  an  amplitude  which 

concordance 

is  still  small  compared  to  the  zav  of  (4?).  In  this  case 

i 

the  instabilities  oi  natural  origin  dominate  over  those 
v/hich  have  their  source  in  machining  irregularities.  They 
dominate  not  only  at  t = tconcordance;  but  thereafter. 
Consequently — provided  that  the  statistical  model  is  valid — 
the  imperfections  in  the  fabrication  can  be  neglected  in 
the  analysis  of  Taylor  instability  in  case  (2).  In  this 
case  there  would  seem  to  be  little  point  in  trying  to 
improve  the  perfection  of  manufacturing.  Money  can  be  saved 
by  loosening  up  on  the  machining  tolerances! 

(3)  The  disturbances  which  have  grown  from  initial  irregulari- 
ties of  reduced  wave  lengths  % and 

fastest  growing 

* machining  have  “"^able  amplitudes  at  t = toon<,oraanoe . 
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From  this  circumstance  it  follows  that  the  machining  toler- 
ances have  been  widened  to  the  critical  limits 

( ^ z ^machining  = ( ^ z ^critical 

beyond  which  any  further  decrease  in  the  quality  of  the 

surface  will  shorten  the  time  to  breakthrough.  This  criti- 
cal limit  is  the  fabrication  criterion  sought  from  the 
present  analysis. 

The  concept  of  a critical  magnitude  for  fabricational  irregulari- 
ties, while  easy  to  describe  in  qualitative  terms,  is  difficult  to 
make  precise  in  the  present  state  of  knowledge.  To  illustrate 
this  point,  it  is  sufficient  to  make  one  very  literal  minded  cal- 
culation of  (£  z)crit  along  the  lines  illustrated  in  Fig.  7,  and 
then  to  look  at  the  uncertainties  that  are  associated  with  this 
calculation.  In  the  figure,  the  fingers  of  natural  origin  are 
taken  to  increase  as  in  the  provisional  statistical  theory  of 
bubbler  merger. 

z = (kj/2)  gt2 

dz/dt  - k^  gt  (49) 

The  rate  of  rise  given  by  (49)  at  the  time 

*0  - * ^/s)* 

(=47  ( *maoh/s)*  for  1c  - 0.02)  (50) 

has  the  value 

(dz/dt )0  = 0.48  (3.83  %machS ).  (51) 
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Fig.  7.  Competition  between  irregularities  of  fabricational  origin 
and  of  natural  origin  for  control  of  Rayleigh-Taylor  instability. 

Top:  Fabricational  irregularities  dominate.  Middle:  Departures 

from  Ideality  of  natural  origin — governed  by  surface  tension  or 
viscosity — dominate.  Bottom:  Irregularities  of  fabricational 

and  natural  origin  contribute  comparably  to  the  later  stages  of 
instability.  The  conditions  in  this  case  determine  the  critical 
magnitude,  dorm  to  which  it  is  desirable  to  reduce  the  fabricational 
irregularities,  but  below  which  any  further  reduction  gives  a de- 
creased or  even  negligible  payoff.  Though  the  general  ideas  are 
clear,  and  the  curves  in  the  diagram  have  been  calculated  on  the 
basis  of  reasonably  well  defined  assumptions,  the  details  of  these 
assumptions  are  not  determined  at  the  present  time  with  any  precision 
by  either  available  observations  or  existing  theory.  Moreover, 
small  changes  in  these  assumptions  make  a big  change  in  the  order 
of  magnitude  of  the  critical  magnitude  (Sz)  ...  for  the  irregulari- 
ties ( $ z )maGh  of  fabricational  origin  at  th8rstart  of  the  accelera- 
tion. Therefore,  the  particular  assumptions  used  in  computing  the 
curves  in  the  diagram  are  recapitulated  here:  (l)  The  scale 


^■fastest  growing  the  natural  irregularities  at  the  start 
is  negligible  compared  to  the  scale  X mach  of  fabricational 
marks.  (2)  The  natural  instabilities  subsequently  have  a canoni- 
cal distribution  in  sizes,  of  which  the  scale  but  not  the  form 
changes  with  timepdue  to  finger  mergers,  with  an  average  height 
z&Y  = (k-j_  /2)  g given  by  a "breakthrough  coefficient"  k-^ 
which  is  arbitrarily  equated  to  0.02  (3)  The  fabricational 

irregularities  grow  exponentially  following  the  Rayleigh-Taylorj^ 
formula  z = ( S z)moh  exp  (g/  2: mach)'2'  t up  to  z = 0.48(3. 83) 2 *mach 
and  thereafter  (with  continuity  of  z and  dz/dt  appoint  A)  grow 
linearly  at  the  rate  (dz/dt)  = 0.48(3. 83g  /Cmach)2*  11118  simplifying 
assumption  neglects  merger  of  the  fingers  resulting  from  fabrica- 
tion until  these  fingers  have  been  overtaken  by  the  fingers  which 
have  grown  from  instabilities  of  natural  origin. 


i 


The  number  agrees  with  the  rate  of  linear  rise  of  fingers  of  radius 
3.83  X , or  of  reduced  wave  length  X , (point  C in  the 
diagram).  Calculate  back  along  this  linear  curve  of  rise  of  the 
fingers  of  fabricational  origin  to  the  place  (point  A in  diagram) 
where  the  straight  line  agrees  in  magnitude  and  slope  with  the 
Rayleigh-Taylor  exponential  curve: 

(dz/dt)A  = 0.48  (3.83  *maoh  g)  i (52) 

ZA  = (l/*)(dz/dt)A  = ( 5:raaoh/e)5  (dz/dt)  (53) 

= (0.48X3.83 )i*,„ach  = 0.94  * mach; 


tA  = tc  " <V  Zi\ )/(dz/dt )yf fWe>)<' 

« (\ach/s>4  l 1 + (0.48X3.83^/2^]  (54) 

l=SI>-5  ( 3lma0i/s)4  for  kl  “ 0-02  1 


Follow  the  exponential  curve  back  to  the  time  t = o to  find  the 
critical  amplitude  for  fabricational  irregularities: 

<*ZW  ■ (0.48X3.83)*  Xmaohe-^ 

= 0.94  X u e-1-  (0.48) (3.83 )?/2ki  (55) 

'mach 

( “ 0-94  >maoh  e-2lt'5  = 2.2  x 10'n  for  lq.  =0.02) 


If  a value  k^  = 0.04  is  adopted  for  the  breakthrough  coefficient 
(slab  of  fluid  accelerable  to  only  25  times  its  thickness  before 
the  fingers  of  natural  origin  penetrate)  then  the  calculated  value 
of  the  critical  amplitude  is  increased  from  (55)  to 

( <f  z)Crit  = 0.94  /^mach  e 75=10  X maCh  (^=0.04^^ 
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A comparably  drastic  change  in  the  order  of  magnitude  of  the  cal- 
culated ( [ z)  will  be  made  by  allowing — as  Fig. 7 does  not — 
crit 

for  some  preliminary  merger  of  the  fingers  which  grow  out  of  dis- 
turbances of  reduced  wave  length  ^ = 3t  • Evidently  it  is  not 

mach 

possible  from  present  information  to  state  a reliable  figure  for 
the  critical  amplitude  (£  z)crit  or  for  the  corresponding  dimen- 
sionless ratio  ($  z ) crit/  /Cniach*  Experiments  to  test  the  validity 
of  the  concept  of  a critical  amplitude  and  to  measure  its  value 
would  therefore  seem  in  order.  They  would  be  a natural  extension 
of  any  observational  program  to  search  for  a breakthrough  coefficient. 

In  Section  II  electrolysis  was  mentioned  as  a technique  to 
build  up  the  proposed  laminar  structure.  It  may  not  be  out  of 
place  to  think  again  of  electrolysis  as  a general  technique  to 
produce  surfaces  of  very  high  perfection.  Such  perfection — 
according  to  analysis  of  Sections  III  and  IV  of  this  report — will 
not  prevent  Taylor  instability,  and  could,  at  the  worst,  be  a 
form  of  "abacadabra"  which  costs  much  time  and  money;  but  it  could 
also  help  to  postpone  the  breakthrough  of  spikes  and  bubbles. 

Electro -polishing  today  is  so  greatly  developed  as  a technology 
for  the  manufacture  of  the  chrome  parts  of  automobiles  and  for  many 
other  purposes  that  it  is  hardly  necessary  to  recall  the  principle 
of  the  polishing: 

(l)  Alternate  cycles  of  current  with  different  and  pro- 
grammed magnitudes  for  the  product  of  (current  density)  • (time) 
in  the  forward  cycle  and  the  reverse  cycle. 
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(2)  Deposition  of  a very  thin  layer  of  material  in  each 
forward  cycle. 

(3)  Inevitable  formation  of  small  bumps  and  spikes  in  this 
process. 

(4)  Concentration  of  lines  of  force  and  field  gradient 
in  the  electrolytic  fluid  near  these  projections. 

(5)  Cessation  of  the  deposition  before  this  concentration 
effect  has  succeeded  in  magnifying  substantially  these  projections. 

(6)  Reversal  of  the  current  in  the  next  phase  of  the  pro- 
gram to  take  advantage  of  this  concentration  of  current  density 
selectively  to  remove  the  projections. 

(7 ) A new  cycle  of  buildup  (one  sign  of  the  product  of 
(current )• (time)  and  electro-deposition  (smaller  magnitude  and 
opposite  sign  for  the  product  of  (current)  • (time). 

The  technology  of  electro-polishing,  however  fully  developed 
for  making  shiny  fenders,  is  not  known  by  us  ever  to  have  been 
used  for  scientific  purposes  by  anyone  with  the  talents  and  passion 
for  precision  which  Henry  A.  Rowland^)  showed  in  the  design  and 
construction  of  an  almost  perfect  screw  for  his  engine  to  rule 
gratings.  Among  the  factors  that  might  be  considered  if  Rowland's 
model  we re  to  be  followed,  this  time  to  construct  an  almost  perfect 
surface,  are  these: 

(l)  Preparation  of  an  almost  ideal  substrate  upon  which 
to  start  the  electro -deposition.  In  this  connection  one  thinks 
particularly  of  surface  tension  plus  the  pressure  of  a thin 
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filling  gas  as  a mechanism  to  produce  a thin  cylinder  of  fluid 
connecting  two  separated  hemispherical  metal  caps. 

(2)  Design  of  a system  of  electrodes,  or  a system  of 
automatic  control  (scanning  by  light  source  plus  photo  electric 
cell  feeding  into  circuit  which  controls  electrolytic  current)  0^ 
both  so  as  to  maintain  uniformity  of  thickness  as  closely  over 
large  distances  as  over  the  very  small  distances  at  which  electro- 
polishing  is  relevant. 

(3)  Control  of  temperature  and  convection  currents  to  a 
completeness  sufficient  to  prevent  imperfections  from  these 
sources  greater  than  a preassigned  magnitude. 

(4)  Method  of  mounting  foil  and  inserting  it  in  final 
device  which  preserves  its  near  perfection  to  the  moment  of  use. 

Surfaces  prepared  in  this  way  might  be  expected  (a)  to  be  most 
suitable  for  the  experimental  investigation  of  the  statistical 
bubble  theory  of  Section  III,  which  focusses  attention  on  irregu- 
larities at  a microscopic  - not  fabricational  - scale  (b)  to  be 
the  most  suitable  for  those  technical  applications  where  it  is 
important  to  minimize  in  any  way  that  is  feasible  the  effects  of 
Taylor  instability. 
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